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SUMMARY 


The  purpose  o£  this  work  is  to  optimize  the  phase- 
difference  polarization  filter  for  use  in  extracting  short  - 
period  regional  data  from  seismic  noise  and  signal  coda. 

The  phases  to  be  extracted  are  the  bodywaves  Pn,  Pg,  Sn,  and 
Sg,  and  the  surface  waves  Lg  and  Rg.  The  data  available  for 
this  work  were  from  events  with  epicenters  in  northeastern 
China,  and  in  Honshu  and  Hokkaido,  Japan.  The  data  used  were 
recorded  at  the  three-component  short -period  site  of  the 
Korean  Seismic  Research  Station  (KSRS) .  The  polarization  fil¬ 
ter  yields  signal-to-noise  ratio  gains  of  approximately  10  dB 
for  regional  compressional  and  shear  waves  (Pn,  Pg,  Sn,  and 
Sg)  and  of  approximately  8  dB  for  regional  surface  waves  (Lg 
and  Rg) .  Direct  estimates  of  ime  detection  capability  for 
regional  P  and  S  (A=10°  to  13°)  made  before  and  after  polar¬ 
ization  filtering  show  that  the  polarization  filter  yields 
detection  improvements  of  approximately  0.6  m^- units  for 
these  phases.  Finally,  the  polarization  filter  can  be  used 
to  separate  mixed  signals,  to  identify  phases,  and  to  extract 
regional  phases  from  seismic  noise. . -  j 


Neither  the  Advanced  Research  Projects  Agency  nor  the 
Air  Force  Technical  Applications  Center  will  be  responsible 
for  information  contained  herein  which  has  been  supplied  by 
other  organizations  or  contractors,  and  this  document  is 
subject  to  later  revision  as  may  be  necessary.  The  views 
and  conclusions  presented  are  those  of  the  authors  and  should 
not  be  interpreted  as  necessarily  representing  the  official 
policies,  either  expressed  or  implied,  of  the  Advanced  Re¬ 
search  Projects  Agency,  the  Air  Force  Technical  Applications 
Center,  or  the  US  Government. 
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SECTION  I 
INTRODUCTION 


A.  SUMMARY  OF  PREVIOUS  WORK 

This  research  effort  is  intended  to  refine  and  test  a 
technique  to  extract  and  identify  the  short-period  regional 
phases  Pn,  Pg,  Sn,  Sg,  Lg,  and  Rg.  The  technique,  called 
phase-difference  polarization  filtering,  was  originally  de¬ 
veloped  for  long-period  teleseismic  data  (Strauss,  1978). 
However,  with  the  move  toward  utilization  of  three-component 
single-site  stations,  and  with  the  increased  interest  in  re¬ 
gional  seismic  data,  it  becomes  desirable  to  apply  the  tech¬ 
nique  to  short-period  regional  phases. 

An  earlier  form  of  polarization  filtering,  used  to  ex¬ 
tract  compress ional  and  shear  waves,  operated  on  the  expect¬ 
ed  angles  of  incidence  of  these  bodywaves  (Lane,  1977;  Strauss, 
1978).  When  tested  on  long-period  teleseismic  data,  however, 
the  performance  of  this  filter  proved  to  be  regionally  depen¬ 
dent;  it  showed  little  change  in  detection  capability  for 
Kurile -Kamchatka  events  as  recorded  at  Guam  and  0.3- 0.4  m^ 
units  improvement  for  the  same  events  as  recorded  at  Mashhad, 
Iran . 


Practical  considerations  dictate  that  the  first  step  in 
designing  regional  short-period  polarization  filters  should 
be  to  use  only  universally  applicable  criteria.  Such  an  ap¬ 
proach  is  described  below. 
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Specifically,  the  phase-difference  polarization  filter 
developed  for  regional  phases  is  a  logical  extension  of  the 
three -component  adaptive  processor  designed  by  Lane  (1973, 
1976)  .  This  processor  searches  for  long-period  waveforms 
which  display  a  90°  phase  difference  between  the  radial  and 
vertical  components  of  motion,  as  is  exhibited  by  Rayleigh 
waves  (retrograde  elliptical  waveforms).  To  extract  com- 
pressional  waveforms,  the  sought-for  phase  difference  is  set 
to  0°,  while  to  extract  shear  waveforms,  the  sought-for  phase 
difference  is  set  to  180°.  To  extract  prograde  elliptical 
waveforms,  the  sought-for  phase  difference  is  set  to  270°. 

A  separate  algor ithm, which  passes  only  that  energy  which  is 
within  a  specified  range  about  the  predicted  arrival  azimuth, 
is  used  to  extract  Lg  data  recorded  on  the  transverse  compo¬ 
nent.  A  detailed  description  of  this  filter  technique  is 
presented  in  Section  II. 

Previous  applications  of  the  phase-difference  polariza¬ 
tion  filter  to  long-period  teleseismic  data  gave  a  0.5  m^  unit 
improvement  (approximate)  in  detection  capability  for  P,  S, 
and  LR  waveforms  (Strauss,  1978).  Preliminary  testing  of  this 
filter  on  short-period  regional  data  indicated  that  the  same 
improvement  in  detection  capability  might  be  expected.  In 
addition,  preliminary  tests  indicated  that  the  polarization 
filter  may  be  useful  in  waveform  identification  (e.g.,  a  wave¬ 
form  passed  by  the  P  model  of  the  filter  can  be  more  confi¬ 
dently  identified  as  a  P  wave). 
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B.  GOALS 


This  report  presents  the  results  of  research  on  the  ap¬ 
plication  of  the  phase-difference  polarization  filter  to  re¬ 
gional  short-period  data.  The  specific  goals  of  this  study 
were: 

•  To  design  three -component  adaptive  filters  which 
extract  short-period  regional  phase  signals; 

•  To  estimate  the  effects  of  such  filters  on  the  de¬ 
tection  capability  of  the  phases  Pn,  Pg,  Sn,  Sg, 

Lg,  and  Rg; 

•  To  test  the  filters  on  seismic  events  with  epicen- 
tral  distances  between  0°  and  20°  from  the  receiver 
site; 

•  To  outline  future,  feasible  extensions  of  polariza¬ 
tion  filter  use  to  the  problem  of  identifying  sig¬ 
nificantly  refracted  or  wave-guided  phases  which 
are  observed  at  regional  distances  (this  pertains 

to  source  discrimination  by  correlation  with  source- 
region-to-station  synthetics). 

To  accomplish  the  first  goal,  the  existing  long-period 
teleseismic  polarization  filter  was  adapted  to  process  short- 
period  regional  data.  The  intent  is  to  find  those  values  of 
the  polarization  filter  program  parameters  which  produce  maxi¬ 
mum  signal-to-noise  ratio  (SNR)  gains  consistent  with  relative¬ 
ly  low  variance  in  these  gains  from  event  to  event.  This  work 
is  described  in  Section  II. 

Work  performed  toward  achieving  the  second  goal  was 
directed  toward  optimizing  the  performance  of  polarization 
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filters  by  determining  the  improvement  in  signal-to-noise 
ratio  produced  by  application  of  the  polarization  filter. 

This  was  accomplished  by  applying  the  filter  to  a  suite  of 
synthetic  seismograms.  In  this  approach,  a  test  signal  is 
scaled  and  added  to  noise  in  order  to  produce  a  synthetic 
signal  of  known  peak  signal-to-RMS-noise  ratio  (SNR).  Vary¬ 
ing  the  scaling  allows  one  to  create  a  suite  of  signals  which 
differ  only  in  SNR.  Processing  such  a  suite  of  data  with  the 
polarization  filter  then  permits  accurate  gauging  of  SNR  im¬ 
provement  which  can  be  achieved  through  application  of  the 
filter.  More  significantly,  detection  capability  improve¬ 
ment  is  accurately  gauged  by  noting  the  input  SNR  of  the 
smallest  detectable  signal  before  and  after  filter  applica¬ 
tion. 


The  test  data  available  for  this  task  consisted  of  three- 
component  short-period  data  recorded  at  the  Korean  Seismic 
Research  Station  (KSRS) .  The  majority  of  the  regional  events 
for  which  data  were  available  had  epicenters  near  Honshu, 

Japan  (A=10°  to  13°).  A  source  region  in  northeastern  China 
(A~8°)  provided  data  for  approximately  50  events.  Relatively 
little  data  were  available  for  epicentral  distances  less  than 
7°.  Thus,  the  data  available  were  sufficient  to  allow  an 
estimation  of  detection  capability  improvement  produced  by 
application  of  the  polarization  filter  for  P  and  S  phases 
observed  at  regional  distances.  These  direct  estimates  of 
detection  capability  improvement  are  used  to  verify  the  esti¬ 
mates  derived  from  the  improvements  determined  using  synthetic 
seismograms.  Work  performed  on  the  second  goal  is  described 
in  Section  III. 
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The  purpose  of  the  third  goal  is  to  illustrate  the 
capability  of  the  polarization  filter  to  extract  and  identify 
phases  and  to  separate  mixed  signals.  This  work  is  described 
in  Section  IV. 

Finally,  the  conclusions  reached  from  this  work  are  pre¬ 
sented  in  Section  V.  This  section  also  presents  recommenda¬ 
tions  for  future  work  involving  the  phase-difference  polariza¬ 
tion  filter  and  other  feasible  extensions  of  the  polarization 
filter  concept. 

It  should  be  noted  that  all  computer  processing  perform¬ 
ed  during  the  course  of  this  task  was  carried  out  in  a  batch 
processing  environment.  Where  pertinent,  the  differences  be¬ 
tween  this  type  of  processing  environment  and  an  interactive 
environment  will  be  noted. 
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SECTION  II 
METHODOLOGY 


A.  DESCRIPTION  OF  THE  METHOD 

The  polarization  filter  described  in  this  report  depends 
on  the  phase  difference  between  the  radial  and  vertical  sig¬ 
natures  for  extracting  compressional ,  shear,  retrograde  el¬ 
liptical,  and  prograde  elliptical  short-period  regional  sig¬ 
nals.  A  separate  algorithm,  described  below,  is  used  to  ex¬ 
tract  Lg.  Figure  I I - 1  summarizes  the  four  models  available 
in  the  polarization  filter. 

Let  the  radial  and  vertical  Fourier  transform  for  any 
given  time  segment  be  represented  by: 


-j“>0t+<i>RO) 

R(«0 

=  Ar(w) 

-jo)  t+4>vO) 

V(oO 

/ — \ 
3 

> 

II 

e  0  v 

where  u>  is  an  arbitrary  reference  frequency  and  4>dO)  and 
4iy(u)  are  the  radial  and  vertical  instantaneous  phases  with 
respect  to  coo.  Thus,  the  radial  and  vertical  components  at 
any  given  frequency  differ  in  phase  by  [^(uj)  -  4>y(u>)]  .  To 
extract  the  short-period  regional  waveforms,  the  test  statis¬ 
tic 
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is  computed.  If  4>0(io)  is  less  than  some  pre-set  value  aQ,  a 
filter  weight  of  '1'  is  assigned  to  the  data.  If  4>0(w)  is 
greater  than  ao  but  less  than  some  value  a^,  a  filter  weight 
between  zero  and  '1'  is  assigned;  the  value  of  the  filter 
weight  depends  on  the  separation  between  <J>oCw)  and  aQ.  If 
<t>o(<j)  is  greater  than  a  filter  weight  of  zero  is  assigned. 
Thus,  the  pass  window  is  trapezoidal  in  shape  (see  Figure 
II-2).  This  procedure  allows  the  acceptance  of  signals  whose 
phase  difference  does  not  exactly  meet  the  appropriate  model 
because  of  noise-induced  distortions.  A  pictorial  overview 
is  presented  in  Figure  1 1  - 3 ;  a  detailed  description  of  the 
polarization  filter  program  is  presented  in  Appendix  A. 

Note  that  the  amount  of  noise  suppression  to  be  realized 
can  be  predicted  for  this  filter.  For  noise  with  random 
radial -vertical  phase  differences,  <J>0  (w)  will  range  between 
0°  and  180°.  Values  of  15°  for  aQ  and  imply  that  one- 

eighth  of  the  noise  will  be  passed  by  the  polarization  filter 
((ao  +  a^/2)/180,  since  the  filter  weights  between  otQ  and 
decrease  linearly  toward  zero).  Rejection  of  seven-eights 
of  the  noise  corresponds  to  a  9  dB  noise  suppression. 

The  concept  of  polarization  filtering  may  further  be 
described  by  considering  signals  and  noise  as  a  superposition 
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FIGURE  II- 3 

PICTORIAL  REPRESENTATION  OF  POLARIZATION  FILTER 
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of  modulated  amplitude  and  phase  components.  The  radial  and 
vertical  time  traces  are  given  as  follows: 


r(t) 

v(t) 


r  -j  W.t  +  4>R(t,0),) 

r  "j  t  +  <J>V(t  ,0).  ) 

•  Ay(t  ,<!)•)  e  1  1 


If  at  any  time  a  regional  phase  occurs  where  the  radial  and 
vertical  component  exceeds  noise  on  one  or  more  of  the  spec¬ 
tral  components,  then  the  desired  phase  relationship  is  sat¬ 
isfied  and  the  signal  component  is  passed  without  distortion. 
The  threshold  for  satisfying  the  phase  relationship  must  be 
set  wide  enough  to  pass  imperfect  signals  and  yet  narrow 
enough  to  reject  chance  correlations  of  random  phase  noise. 

For  example,  given  a  phase  relationship  for  P  waves  of  zero, 
experience  may  indicate  that  at  least  a  ±15°  passband  is  need¬ 
ed  to  pass  signals  without  distortion  Oo  *  =  15°,  in  this 

case).  Further, assume  that  noise  is  equally  likely  to  ex¬ 
hibit  any  phase  angle  between  0°  and  360°.  Then  the  prob¬ 
ability  of  noise  passing  the  signal  criteria  is  1/12.  Thus, 

only  1/I2t^1  of  the  spectral  noise  components  are  weighted  by 

t  s 

one;  the  other  11/12  are  weighted  zero.  By  application 
of  Parseval's  theorem  (Brand,  1955)  the  variance  of  noise  in 
the  time  domain  is  reduced  to  1/I2t^1  of  what  it  was  before  ap¬ 
plying  the  phase  angle  criteria  of  polarization  filters.  There¬ 
fore,  the  theoretical  gain  expected  by  polarization  filters  in 
this  case  is  10  log^  (1/12)=1Q .  8  dB.  An  interesting  result  of 
this  very  simple  noise  model  is  that  this  gain  of  0.5  to  0.6 
m^  units  of  improvement  would  be  expected  for  any  propagation 
mode  to  which  this  criteria  is  applied;  the  reason  for  this 
is  that  the  random  phase  of  the  noise  is  independent  of  the 
propagation  mode  to  be  extracted.  Experience  (Strauss,  1977) 
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with  long-period  surface  wave  extraction  as  well  as  short- 
period  regional  phase  extraction  indicates  this  to  be  the 
case.  Some  variation  from  this  expected  gain  might  be  en¬ 
countered  in  cases  where  the  radial  and  vertical  components 
differ  greatly.  Fot  example,  if  Ay>>AR  as  for  teleseismic  P 
waves,  higher  signal-to-noise  ratios  would  be  required  for 
both  components  to  dominate  the  noise.  In  such  a  case,  the 
expected  gain  would  be  considerably  less  than  predicted.  For 
long-period  surface  waves  and  the  regional  phases  extracted 
by  this  criteria,  Ay~AR. 


To  extract  Lg,  the  assumption  is  made  that  Lg  is  the 
short-period  analog  of  the  long-period  Love  wave.  As  shown 
in  Figure  I I- 4,  the  apparent  arrival  azimuth  of  the  Love  (Lg) 
wave  can  be  expected  to  show  small  deviations  (3)  from  the 
arrival  azimuth  predicted  from  the  event  epicenter  and  re¬ 
ceiver  location.  Noise,  however,  will  produce  values  of  3 
varying  randomly  from  0°  to  360°.  The  algorithm  used  to 
extract  Lg  waves,  therefore,  computes  this  angle  3  on  a 
frequency-by-frequency  basis  using 

-1  Rr 

3  =  tan  sri-  , 
xr 

where  Rr  is  the  real  part  of  the  Fourier  transform  of  the 
radial  component  while  T  is  the  real  part  of  the  transform 
of  the  transverse  component. 


The  computed  value  of  3  is  used  to  assign  a  filter 
weight  to  the  Fourier  amplitude  at  the  frequency  from  which 
it  was  computed,  where  the  filter  weight  is  determined  in 
the  same  fashion  as  previously  described.  This  algorithm  is 
discussed  in  detail  by  Strauss  (1976). 
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FIGURE  1 1 - 4 

LOVE  WAVE  ARRIVAL  AZIMUTH  FOR  AN 
ARBITRARY  FREQUENCY  COMPONENT 
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B.  POLARIZATION  FILTER  OPTIMIZATION 


It  now  remains  to  optimize  the  performance  of  the  polar¬ 
ization  filter  (in  terms  of  producing  the  greatest  SNR). 
Optimization  of  the  filter  will  free  the  analyst  from  the 
necessity  to  repeatedly  process  each  event  in  order  to  ob¬ 
tain  the  'best'  result.  In  optimizing  the  use  of  the  polar¬ 
ization  filter,  minimal  prior  knowledge  of  the  signal  to  be 
processed  is  assumed  (all  that  must  be  known  about  the  data 
is  the  source  location  since  the  filter  does  require  data 
rotated  to  a  vertical,  transverse,  radial  configuration). 

It  has  been  demonstrated  (Lane,  1977)  for  long-period 
signal  extraction  applications  that  a  bandpass  prefilter 
produces  greater  SNR  improvements  than  does  application  of 
a  polarization-type  filter  alone.  Whether  this  holds  true 
in  the  case  of  short-period  regional  data  must  be  examined, 
and,  if  it  does,  the  question  of  what  bandpass  filter  serves 
as  the  best  prefilter  must  be  answered. 

Table  1 1 - 1  shows  the  results  of  processing  a  suite  of 
regional  P  waves  and  a  suite  of  regional  S  waves  (A~10°)  us¬ 
ing  the  polarization  filter.  Each  data  set  was  processed 
once  with  the  polarization  filter  applied  to  the  raw  data, 
and  once  with  the  polarization  filter  applied  to  bandpass- 
filtered  data.  As  is  seen,  the  bandpass/polarization  filter 
cascaded  combination  produces  higher  SNR  gains  with  lower 
variance  in  the  gains  than  does  the  polarization  filter  alone. 
Thus,  it  appears  desirable  to  preprocess  data  using  a  bandpass 
filter  prior  to  applying  the  polarization  filter. 

The  bandpass  filter  used  in  the  above  test  was  a  broad¬ 
band  type  having  a  passband  of  0. 5-5.0  Hz  and  a  10-point 
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cosine-squared  taper  (data  sampled  at  10  samples  per  second) . 

It  is  always  possible  to  achieve  higher  SNR  gains  by  tailor¬ 
ing  the  passband  to  the  particular  waveform  to  be  extracted. 

For  example,  if  it  is  known  that  Sn  from  a  particular  source 
region  is  dominated  by  2  Hz  energy  at  a  given  receiver,  a 
narrowband  filter  centered  on  this  frequency  will  produce 
higher  SNR  gains  than  will  a  0. 5-5.0  Hz  bandpass  filter. 

There  is,  however,  a  danger  in  tailoring  the  filter  to  the 
expected  dominant  frequency  of  the  waveform  since  this  re¬ 
quires  prior  knowledge  of  the  spectrum  of  each  phase  from 
each  region  at  each  receiver.  To  see  this,  consider  Figure 
1 1 -  5 ,  which  shows  the  shear  wave  (Sn-Sg)  gate  of  an  event 
from  northeastern  China  which  has  been  processed  by  the  polar¬ 
ization  filter,  but  which  has  been  preceded  by  three  different 
bandpass  filters.  Each  pair  of  traces  in  this  figure  consists 
of  the  vertical  trace  after  application  of  the  bandpass  filter 
(upper  trace  of  each  pair),  and  the  same  trace  after  applica¬ 
tion  of  the  bandpass  filter  and  the  shear  model  of  the  polar¬ 
ization  filter.  The  scales  (S)  given  are  the  maximum  ampli¬ 
tude  of  each  trace  in  digital  counts.  Also  shown  are  the 
SNRs  for  Sn  and  Sg,  as  appropriate.  This  figure  indicates 
that  use  of  an  appropriately-tailored  narrowband  filter  prior 
to  polarization  filtering  yields  higher  SNR  gains  than  is 
possible  using  a  broadband  preprocessing  filter.  However, 
note  that  application  of  the  1.25-5.0  Hz  narrowband  filter 
(center  pair  of  traces)  eliminates  Sg  while  it  enhances  Sn . 
Further,  application  of  the  0.5-1.25  Hz  narrowband  filter 
(bottom  pair  of  traces)  eliminates  Sn  while  it  enhances  Sg. 

In  comparison  to  this,  the  wideband  filter  (top  pair  of  traces) 
enhances  (albeit  to  a  lesser  degree)  both  Sn  and  Sg .  Thus, 
when  the  polarization  filter  program  is  used  in  a  batch  pro¬ 
cessing  mode,  it  appears  preferable  to  use  a  wideband 
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FIGURE  1 1  -  5 

COMPARISON  OF  THREE  BANDPASS  FILTERS  CASCADED  WITH  THE 
PHASE- DIFFERENCE  POLARIZATION  FILTER  (VERTICAL  COMPONENT) 


preprocessing  filter.  For  the  remainder  of  the  work  report¬ 
ed  here  (unless  specifically  noted),  the  bandpass  filter 
used  has  a  0. 5-5.0  Hz  passband  with  a  10-point  cosine-squared 
taper. 

It  should  be  noted  that  the  use  of  several  narrowband 
filters  would  be  of  great  value  in  an  interactive  processing 
environment.  Here,  the  analyst  could  select  separate  time 
gates  for  each  phase.  He  could  then  make  repeated  applica¬ 
tions  of  the  polarization  filter,  using  a  different  bandpass 
filter  for  each  run  in  order  to  determine  the  best  passband 
to  cascade  with  the  polarization  filter  for  each  phase.  As 
the  work  progressed,  the  analyst  would  learn  the  best  pass- 
band  for  each  phase,  and  would  know  whether  source  region  or 
source-receiver  path  effects  impact  on  the  choice  of  a  pass- 
band. 


The  next  point  to  be  considered  is  the  question  of  how 
the  controlling  parameters  of  the  polarization  filter  are  to 
be  set . 

The  following  parameters  control  the  operation  of  the 
polarization  filter  algorithm:  the  number  of  points  per 
processing  segment,  and  the  two  parameters  which  determine 
the  width  of  the  trapezoidal  pass  window  (aQ  and  ,  as  de¬ 
scribed  earlier)  .  Since  filter  weights  are  computed  and  ap¬ 
plied  on  a  frequency-by-frequency  basis  for  each  segment,  in¬ 
creasing  the  segment  length  for  a  finite  signal  length  has 
the  effect  of  increasing  the  effect  of  noise;  that  is,  it 
causes  distortion  of  the  signal  phase  difference.  At  some 
point,  the  segment  length  will  be  so  long  that  the  signal  in 
the  segment  will  be  rejected  along  with  the  noise.  On  the 
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other  hand,  too  short  a  segment  length  results  in  a  less- 
than-optimum  result.  To  see  this,  consider  the  case  where 
the  first  few  points  of  a  signal  lie  near  the  end  of  a  data 
segment.  The  greater  part  of  the  segment  is  noise,  and  so, 
these  first  few  points  of  the  signal  will  be  rejected  by  the 
filter  as  noise.  Although  the  next  segment  will  contain  the 
signal  (normally  a  50-percent  segment  overlap  is  used),  the 
initial  portion  of  the  signal  will  be  distorted  due  to  the 
earlier  rejection  of  the  beginning  of  the  signal.  Thus,  it 
is  necessary  to  estimate  the  segment  length  at  which  this 
problem  is  minimized. 

The  effects  of  varying  the  width  of  the  phase  pass  win¬ 
dow  are  somewhat  easier  to  state.  The  radial -vertical  phase 
difference  of  any  signal  will,  to  some  extent,  always  be  dis 
torted  by  the  underlying  noise.  Widening  the  pass  window 
allows  more  of  the  signal  to  pass;  however,  such  widening 
also  allows  more  noise  to  pass.  In  any  given  case,  there 
should  be  a  width  such  that  the  signal  is  passed  with  a  mini 
mum  of  accompanying  noise. 

The  approach  taken  to  estimate  the  best  values  for  the 
pass  windows  for  the  various  phases  (in  terms  of  maximizing 
the  SNR)  was  empirical.  Specifically,  a  suite  of  test  data 
for  each  phase  was  scaled  and  added  to  noise  recorded  just 
prior  to  the  first  arrival  of  the  signals.  The  resulting 
traces  had  SNRs  so  low  that  the  signals  could  barely  be  seen 
in  the  noise  background.  The  test  data  were  then  repeatedly 
processed  using  the  polarization  filter,  with  different 
values  of  the  program  parameters  used  for  each  run.  The 
test  values  selected  for  the  number  of  points  per  processing 
segment  were  16,  32,  and  64,  while  the  test  values  initially 
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selected  for  the  pass  window  half-width  parameters  (aQ=a^) 
were  5°,  10°,  15°,  and  20°.  Note  that  preliminary  tests 
revealed  that  the  pass  window  half -width  of  5°  produced  very 
large  variances  in  the  mean  signal-to-noise  ratio  values; 
thus,  this  test  value  for  the  pass  window  half-width  was 
dropped  from  further  consideration. 

The  results  of  this  test  are  summarized  in  Table  II-2. 

Note  that  in  each  test,  the  data  were  prefiltered  with  a 
0. 5-5.0  Hz  bandpass  filter  before  application  of  the  polar¬ 
ization  filter.  The  phase  Lg  was  not  separately  tested  since 
the  polarization  program  operates  on  Lg  and  Rg  in  the  same 
pass.  Thus,  the  same  program  parameters  are  used  for  both 
phases.  The  phase  Pg  was  not  tested  since  insufficient  data 
were  available  to  make  this  test.  However,  in  the  batch  pro¬ 
cessing  environment  used  throughout  this  study,  the  parameters 
used  to  process  Pn  were  also  _sed  to  process  Pg.  The  data  of 
Table  1 1 - 2  show  that  16  points  per  processing  segment  yielded 
higher  SNR  gains  than  did  segment  lengths  of  32  or  64  points. 

It  is  instructive  to  examine  pictorially  the  effect  on  the 
polarization  filter  output  produced  by  varying  the  number  of 
points  per  processing  segment.  With  reference  to  Figure  II-6, 
the  top  set  of  traces  shows  an  Lg-Rg  time  gate  after  bandpass 
filtering.  The  three  sets  of  traces  below  this  first  set 
show  the  polarization  filter  output  using  16,  32,  and  64 
points  per  processing  segment,  respectively  (the  pass  window 
half-width  was  set  to  15°  in  each  case  (see  below)).  As  the 
processing  segment  length  increases,  the  apparent  length  of 
the  passed  waveform  increases,  thereby  making  more  difficult 
the  picking  of  the  phase  start  time  (i.e.,  its  first  motion). 
Also,  the  peak  amplitude  of  the  waveform  passed  decreases  as 
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SIGNAL- TO -NOISE  RATIO  GAINS  PRODUCED  BY  VARIOUS  COMBINATIONS 
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FIGURE  II-6 

EFFECTS  OF  CHANGING  NUMBER  OF  POINTS  PER  PROCESSING 
SEGMENT  ON  POLARIZATION  FILTER  RESULTS 
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the  processing  segment  increases,  in  most  cases,  indicating 
degradation  of  the  desired  signal  by  interference  with  noise 
and  the  other  waveforms  in  the  Lg-Rg  window.  This  result, 
and  that  described  above,  suggest  that  a  segment  length  of 
16  points  is  preferred. 

The  choice  of  the  width  of  the  phase  pass  window,  how¬ 
ever,  is  not  clear  from  the  results  of  the  data  in  Table  1 1 - 2 . 
Though  a  10°  window  width  appears  to  perform  better  than  the 
other  window  widths  examined,  the  standard  deviation  associated 
with  the  mean  gain  values  suggest  that  the  mean  gain  determina¬ 
tions  are  not  statistically  different.  Therefore,  the  middle 
value  of  the  three  window  widths,  15°,  was  arbitrarily  select¬ 
ed  for  the  phase  pass  window  half-wdith  to  be  used  in  this 
work . 


In  sum,  a  16-point  time  window  and  a  phase  window  half- 
width  of  ao=15°  are  preferred,  and  these  values  were  used  for 
all  processing  described  in  the  remainder  of  this  report. 

The  next  section  examines  the  effects  produced  by  the  Polar¬ 
is  ization  filter  on  SNR  and  detection  capability. 
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SECTION  III 


THE  EFFECTS  OF  POLARIZATION  FILTERING  ON  SIGNAL-TO-NOISE 
RATIO  AND  DETECTION  CAPABILITY 


This  section  examines  the  effect  of  polarization  filter¬ 
ing  on  the  signal-to-noise  ratios  (SNR)  of  the  short-period 
regional  phases  and  the  resulting  improvement  in  detection 
capability  of  these  phases. 

The  'signal  buried  in  noise'  technique  is  used  to  esti¬ 
mate  SNR  improvements  resulting  from  polarization  filtering. 
To  use  this  technique,  one  selects  a  waveform  which  has  a 
high  SNR  and  which  is  preceded  by  seismic  noise.  The  signal- 
centered  portion  of  the  record  is  scaled  by  a  factor  between 
zero  and  '1',  and  is  added  to  the  preceding  noise  (this  pro¬ 
cess  enables  the  user  to  control  the  SNR  of  data  used  to  test 
a  given  signal  extraction  technique).  By  varying  the  scale 
factor,  a  suite  of  three-component  composite  records  can  be 
created  covering  any  desired  range  in  SNR. 

The  first  step  toward  establishing  the  SNR  gain  curve 
for  the  polarization  filter  is  to  apply  the  polarization  fil¬ 
ter  to  each  three-component  composite  record  set.  The  SNR  of 
each  component  is  measured  before  and  after  the  filtering 
process  so  that  the  gain  due  to  polarization  filtering  (i.e., 
the  difference  between  input  and  output  SNR)  can  be  deter¬ 
mined.  In  computing  these  SNR  values,  the  point  at  which  the 
signal  amplitude  is  measured  is  set  to  the  time  of  the  maxi¬ 
mum  value  of  the  original  signal.  This  produces  SNR  measure¬ 
ments  which  are  comparable  even  at  low  input  SNR  values. 
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To  provide  an  objective  measure  of  input  SNR  against 
which  to  plot  the  polarization  filter  gain,  a  quantity  term¬ 
ed  'true  input  signal -to-noise  ratio'  is  computed.  This 
quantity  is  defined  as  follows: 

TRUE  INPUT  SNR  =  20  *  L°g10 

where 

'SF'  is  the  scale  factor,  and 

'Max.  Amp.'  is  the  maximum  zero-to-peak  signal  amplitude 

Figure  1 1 1  - 1  shows  two  examples  of  composite  traces  form¬ 
ed  in  the  manner  described  above  as  they  appear  before  and 
after  polarization  filtering.  The  data  were  created  from  a 
P  wave  recorded  for  an  event  with  an  epicenter  located  8° 
from  the  recording  site.  The  true  input  SNR  of  Set  A  is  ap¬ 
proximately  7  dB  while  the  true  input  SNR  of  Set  B  is  approxi¬ 
mately  15  dB.  The  apparent  gain  in  SNR  produced  by  polariza¬ 
tion  filtering  was  approximately  8  dB  for  the  data  of  Set  A, 
and  was  approximately  12  dB  for  Set  B. 

Application  of  the  polarization  filter  to  a  suite  of 
data  having  true  input  SNR  ranging  from  -4  to  20  dB  yields 
the  example  of  a  typical  SNR  gain  curve  shown  in  Figure  III -2. 
In  this  figure,  the  polarization  filter  gain  appears  to  be 
non-linear  for  true  input  SNR  values  below  approximately  14 
dB.  Above  this  point,  the  gain  curve  flattens,  indicating 
a  constant  polarization  filter  gain  of  12  dB.  This  saturation 
level  gain  has,  in  the  past  (e.g.,  Lane,  1973),  been  taken  as 
the  filter  gain.  The  point  to  be  resolved  is  whether  the  gain 
is  non-linear  below  this  point. 


(SF)  *  (Max.  Amp.) 
RMS  Noise 
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POLARIZATION  FILTER  GAIN  AS  A  FUNCTION  OF  TRUE  INPUT  SIGNAL-TO-NOISE  RATIO 


This  matter  of  apparent  gain  non-linearity  can  be  re¬ 
solved  by  determining  the  difference  between  the  true  input 
SNR  at  which  the  signal  can  be  detected  before  polarization 
filtering  and  the  true  input  SNR  at  which  the  signal  can  be 
detected  after  polarization  filtering.  To  provide  an  objec¬ 
tive  estimate  of  these  values  of  true  input  SNR,  an  automatic 
detector  of  the  type  formulated  by  Unger  (1978)  was  added  to 
the  polarization  filter  program.  The  results  of  this  detec¬ 
tor  showed  that  the  gain  curve  approaches  zero  slope  at  about 
the  lowest  true  input  SNR  (approximately  15  dB)  at  which  the 
signal  can  be  detected  both  before  and  after  polarization 
filtering.  This  point  is  marked  by  ’ B  *  in  Figure  1 1 1  -  2 .  The 
lowest  true  input  SNR  at  which  the  signal  can  be  detected  by 
this  automatic  detector  after  polarization  filtering,  marked 
by  'A'  in  Figure  1 1 1  -  2  ,  is  approximately  4  dB .  These  results 
yield  a  polarization  filter  gain  of  11  dB,  which  agrees  quite 
well  with  the  previously  mentioned  value  of  12  dB . 

This  result  suggests  that  the  polarization  filter  gain 
is  constant.  An  explanation  for  the  apparent  non-linear 
behavior  indicated  by  the  data  of  Figure  III -2,  however,  is 
that  below  the  lowest  true  input  SNR  at  which  the  signal  is 
detected  before  polarization  filtering  (Point  B)  ,  the  gain 
determinations  are  increasingly  influenced  by  noise.  Thus, 
gain  appears  to  decrease  with  decreasing  true  input  SNR. 

This  effect  is  an  aberration  of  the  method  used  to  gauge 
processing  gains  and,  apparently,  it  does  not  accurately  re¬ 
flect  the  performance  of  the  polarization  filter. 

Figure  1 1 1  -  5  shows  the  polarization  filter  gain  versus 
true  input  SNR  curves  for  the  Pn ,  Pg ,  and  regional  (A=12°) 

P  phases.  The  curve  for  Pn  was  determined  from  a  northeastern 
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True  Input  Signal-To-Noise  Ratio  (dB) 

FIGURE  1 1 1- 3 

COMPRESS  I ONAL  WAVEFORM  POLARIZATION  FILTER  GAIN  CURVES  FOR  Pg,  Pn  ,  AND  REGIONAL 


China  event,  while  the  curve  for  Pg  was  determined  from  a 
Kyushu,  Japan,  event;  the  curve  for  regional  P  was  determined 
from  a  Honshu,  Japan,  event.  The  key  data  points  in  this 
figure  are  again  marked  'A'  and  'B';  as  before,  Point  A  marks 
the  lowest  input  SNR  for  which  the  automatic  detector  de¬ 
clared  a  detection  on  the  polarization  filtered  data  while 
Point  B  marks  the  lowest  input  SNR  for  which  the  automatic 
detector  declared  a  detection  on  the  bandpass  filtered  data. 
These  data  (i.e.,  Points  A  and  B)  suggest: 

•  Application  of  the  polarization  filter  improves 
the  SNR  of  compress ional  waves  by  approximately 
10  dB;  this  implies  a  lowering  of  the  compres- 
sional  wave  detection  threshold  of  0.5  m^  units; 

•  Application  of  the  polarization  filter  renders 
detectable  a  compress ional  waveform  with  an  SNR 
as  recorded  (i.e.,  'true')  as  low  as  2-7  dB . 


Figure  1 1 1  -  4  shows  the  polarization  filter  gain  versus 
true  input  SNR  for  Sn  and  regional  (A=12°)  S  phases  (polar¬ 
ization  filter  analysis  of  regional  data  for  the  phase  Sg 
presents  a  sp:cial  problem  which  will  be  discussed  later  in 
this  section) .  The  curve  for  Sn  was  determined  from  a  north¬ 
eastern  China  event,  while  the  curve  for  regional  S  was  de¬ 
termined  from  a  Honshu,  Japan,  event.  The  Points  A  and  B  on 
this  figure  are  as  previously  defined  for  Figure  1 1 1  -  3 . 

These  data  indicate: 


•  The  polarization  filter  increases  the  SNR  of  shear 
waves  (S  or  Sn)  by  approximately  10  dB;  this  im¬ 
plies  a  lowering  of  the  shear  wave  detection  thresh¬ 


old  of  0.5  m^  units; 
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•  Application  of  the  polarization  filter  renders  de¬ 
tectable  a  shear  waveform  with  SNR,  as  recorded, 
as  low  as  2-4  dB. 

A  problem  arises  when  one  attempts  to  determine  the  ef¬ 
fect  the  polarization  filter  has  on  the  phases  Lg ,  Rg,  and 
Sg;  this  is  especially  true  for  cases  where  the  event  epi¬ 
center  is  at  a  distance  of  less  than  10°  from  the  recording 
site.  With  estimated  velocities  of  3.5  km/sec  for  Lg  (Ewing, 
Jardetzkv,  and  Press,  1957),  3.37  km/sec  for  Sg  (Jeffreys 
and  Bullen,  1967),  and  3.2  km/sec  for  Rg ,  one  observes  these 
phases  to  overlap  in  time  at  small  epicentral  distances  (see 
Figure  1 1 1 -  5 ) .  This  overlap  makes  it  difficult  to  identify 
the  maximum  amplitude  of  any  one  of  these  phases  on  the  band¬ 
pass  filtered  traces.  Further,  it  is  not  desirable  to  use 
the  amplitude  observed  at  the  predicted  phase  arrival  time 
since  the  above-mentioned  velocities  are  not  the  same  every¬ 
where  (e.g.,  Lg  velocities  as  low  as  3.2  km/sec  have  been 
noted  by  Herrin  and  Richmond  (I960)).  Given  the  above,  esti¬ 
mates  of  input  SNR  for  Lg ,  Rg,  and  Sg  will  not  be  reliable, 
and  so,  no  curves  of  polarization  filter  gain  versus  input 
SNR  are  presented  for  these  phases. 

Since  it  was  not  possible  to  examine  in  detail  the  ef¬ 
fect  of  the  polarization  filter  on  SNRs  for  Lg,  Rg,  and  Sg , 
as  was  done  for  the  compress ional  phases  and  for  the  shear 
phases  S  and  Sn,  the  polarization  filter  gain  for  these  phases 
were  estimated  in  another  manner.  In  particular,  a  suite  of 
test  events  showing  the  overlapping  Lg,  Sg,  and  Rg  waveforms 
similar  to  those  shown  in  Figure  1 1 1 -  5  was  processed  using 
the  shear  and  surface  wave  models  of  the  polarization  filter. 
Here,  however,  the  effects  of  the  filter  on  noise  and  on  the 
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signal  (as  determined  by  the  estimated  time  of  arrival)  were 
determined  separately.  The  difference  in  the  polarization 
filter's  performance  on  these  waveforms,  then,  is  a  measure 
of  processor  gain.  Suites  of  compressional  and  shear  data 
were  also  processed  in  the  same  manner  so  as  to  provide  gain 
estimates  which  could  be  compared  to  those  previously  deter¬ 
mined. 

The  results  of  this  analysis  effort  are  summarized  in 
Table  1 1 1  - 1 .  The  table  gives  the  sample  population  means  in 
dB  (plus -or -minus  one  standard  deviation)  for  the  degree  of 
noise  suppression,  the  degree  of  signal  suppression,  and  the 
estimated  gain  in  SNR  produced  by  polarization  filtering. 

The  polarization  filter  gains  for  Sg,  Lg,  and  Rg  are  13  dB, 

8  dB,  and  8  dB,  respectively.  It  is  of  interest  to  note  that 
gains  in  SNR  produced  by  the  polarization  filter  for  the  com¬ 
pressional  and  shear  (S  and  Sn)  phases  as  determined  by  this 

4 

method  agree  well  with  the  gains  determined  earlier  by  the 
signal-buried-in-noise  method  (i.e.,  gains  ~10  dB) ,  thereby 
lending  credence  to  the  method  used  to  measure  the  processor 
gains  determined  for  Sg,  Lg,  and  Rg. 

The  SNR  gains  of  Table  III-l  imply  increases  in  the  de¬ 
tection  threshold  of  0.5  m^  units  for  compressional  and  shear 
(S,  Sn)  waves,  0.6  m^  units  for  Sg,  and  0.4  m^  units  for  Lg 
and  Rg. 

Confirmation  of  the  detection  threshold  results  is  sought 
by  making  direct  estimates  of  the  detection  thresholds  before 
and  after  application  of  the  polarization  filter.  The  data 
available  for  this  work,  however,  only  allowed  computation  of 
direct  estimates  of  the  detection  thresholds  for  regional  P 
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TABLE  II I -  1 

EFFECTS  OF  POLARIZATION  FILTER  ON  THE 
SIGNAL -TO -NOISE  RATIOS  OF  REGIONAL  PHASES 


Phase 

Noise 

Suppression 

(dB) 

Signal 

Suppression 

(dB) 

Net  Gain 
(dB) 

Sg 

16.8  +  3 .  S 

4.0  +  3.6 

12.8  +  5.4 

Lg 

10.4  +  1.8 

2.6  +  1.8 

7.9  +  2.2 

Rg 

12.1  +  2.7 

4.0  +  1.6 

8.1  +  3.0 

P,  Pn 

14.0  +  3.2 

3.6  +  1.8 

10.4  +  3.4 

S ,  Sn 

16.6  +  4.3 

6.0  +  4.4 

10.5  +  6.0 

Values  shown  are  sample  population  means  in  dB 
(plus -or-minus  one  standard  deviation). 
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and  regional  S;  insufficient  data  were  available  to  make 
these  estimates  for  Pn,  Pg,  Sg,  and  Sg.  Of  the  approximately 
340  events  studied,  only  about  40  had  epicentral  distances 
such  that  these  phases  might  have  been  detected.  However,  a 
group  of  123  P  waves  and  120  S  waves  from  events  in  the  Hon¬ 
shu  and  Hokkaido,  Japan,  areas  (epicentral  distances  =  10°  to 
13°)  were  available  for  direct  estimates  of  detection  capabil¬ 
ity  for  regional  P  and  S.  Unfortunately,  these  could  not  also 
be  used  for  Lg  and  Rg  since  the  Sea  of  Japan,  which  lies  be¬ 
tween  the  source  and  the  receiver,  blocks  transmission  of  Lg 
(and  perhaps  Rg,  as  well,  as  no  Rg  was  observed  on  any  of 
these  records).  This  situation  is  to  be  expected  since  Ewing, 
Jardetzky,  and  Press  (1957)  note  that  as  little  as  2°  of  ocean 
crust  between  source  and  receiver  will  block  transmission  of 
Lg.  Thus,  it  is  possible  to  make  direct  estimates  of  detec¬ 
tion  capability  thresholds  only  for  regional  P  and  S.  These 
should,  however,  suffice  to  confirm  or  deny  the  previously 
stated  changes  in  detection  threshold  determined  from  the 
computed  signal -to -noise  ratio  gains. 

The  method  used  for  making  direct  estimates  of  detection 
capability  is  to  examine  visually  the  data  traces  recorded  for 
each  event  after  bandpass  filtering  and  again  after  polariza¬ 
tion  filtering  in  order  to  determine  whether  the  waveform  is 
detectable.  Statistics  of  the  detection  status  for  the  events 
of  the  data  base  were  then  compiled  and  were  used  to  compute 
the  percentages  of  events  detected  at  each  bodywave  magnitude 
increment  (bodywave  magnitudes  are  those  reported  by  the 
National  Earthquake  Information  Service).  The  cumulative 
Gaussian  probability  function  is  then  fitted  to  these  detec¬ 
tion  percentage  statistics  by  the  maximum  likelihood  method 
formulated  by  Ringdal  (1974). 
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Detection  threshold  results  are  presented  in  Figures 
1 1 1 -  6  (for  P  waves)  and  1 1 1  -  7  (for  S  waves).  Figures  III-6A 
and  III-7A  show  the  results  for  the  data  after  bandpass  fil¬ 
tering  but  before  application  of  the  polarization  filter; 
Figures  III-6B  and  III-7B  show  the  results  for  the  data  after 
polarization  filtering.  The  upper  portion  of  each  figure 
shows  the  number  of  detections  and  non-detections  at  each 
bodywave  magnitude  increment;  the  lower  portion  shows  the 
resulting  detection  percentages  and  the  fitted  Gaussian 
probability  function.  The  statistic  on  which  the  before  and 
after  polarization  filtering  comparison  is  made  is  the  50 
percent  detection  threshold,  denoted  by  'm^SO'  in  these  fig¬ 
ures.  Comparison  of  the  m^50  values  shows  a  lowering  of  the 
detection  threshold  of  0.59  m^  units  for  P  waves  and  0.5c  m^ 
units  for  S  waves.  These  values  compare  well  with  the  detec¬ 
tion  threshold  improvements  determined  using  the  'signal- 
bur  ied-  in-noise  '  method  (0.5  m^  units  in  both  cases). 

In  summation,  the  data  indicate  that  the  polarization 
filter  may  be  expected  to  improve  the  detection  capability 
of  a  three-component  short-period  station  by  approximately 
0.5  to  0.6  m^  units  for  Pn,  Pg,  and  regional  P,  by  approxi¬ 
mately  0.5  to  0.6  m^  units  for  Sn ,  Sg,  and  regional  S,  and 
by  approximately  0.4  m^  units  for  Lg  and  Rg. 
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SECTION  IV 

EXAMPLES  OF  POLARIZATION  FILTER  USAGE 


This  section  presents  examples  of  polarization  filtered 
data.  These  examples  were  selected  to  demonstrate  the  value 
of  this  filter  in  the  areas  of  signal  extraction,  signal 
identification,  and  mixed  signal  separation.  In  each  figure 
in  this  section,  the  notations  to  the  left  of  each  set  of 
traces  denote  the  following: 

•  BPF  -the  bandpass -filtered  data  which  were  input 

to  the  polarization  filter  program  (0. 5-5.0 
Hz  passband  unless  otherwise  specified) 

•  P  -  the  data  output  by  the  compress icnal  model 

of  the  polarization  filter 

•  Lg,Rg  -  the  data  output  by  the  retrograde  elliptical 

model  of  the  polarization  filter  and  the 
transverse  component  Lg  algorithm 

•  S  the  data  output  by  the  shear  model  of  the 

polarization  filter. 

The  letters  V,  T,  R  indicate  the  vertical,  transverse,  and 

radial  components  of  motion,  respectively.  The  notation 

’S=’(some  number)  indicates  the  maximum  amplitude  of  the 

trace  in  digital  counts.  The  tick  marks  labeled  P,  Pn,  Pg, 

S,  Sn,  Sg,  Lg,  or  Rg  indicate  the  predicted  arrival  time  of 

the  designated  phase.  When  a  symbol  of  the  type  Pn  appears 

'i' 

above  a  trace,  it  indicates  a  visually  picked  phase  arrival 
of  the  type  named  in  the  symbol. 
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Figure  IV-1  illustrates  the  extraction  of  Pn  and  Pg  from 
a  record  with  low  SNR.  Figure  IV-1A  shows  the  results  of  ap¬ 
plying  the  compress ional  model  of  the  polarization  filter  to 
data  which  were  preprocessed  with  the  standard  0. 5-5.0  Hz 
bandpass  filter.  In  the  bandpass -filtered  data,  it  is  barely 
possible  to  detect  the  Pn  arrival,  while  the  Pg  arrival  is 
not  detectable.  After  polarization  filtering,  the  Pn  and  Pg 
phases  are  visually  detectable,  as  indicated.  Figure  IV-1B 
is  presented  to  show  the  improvement  which  can  be  achieved 
by  using  narrowband  bandpass  filter  (1.25-5.0  Hz)  as  a  pre¬ 
processor  to  the  polarization  filter  when  the  dominant  fre¬ 
quency  of  the  sought-after  phases  can  be  estimated.  These 
results  serve  to  demonstrate  that  interactive  processing,  with 
its  ease  of  selection  of  time  gates  and  program  parameters, 
is  more  suitable  for  polarization  filter  application  than  is 
the  batch  processing  environment. 

A  second  example  of  regional  waveform  extraction  is 
shown  in  Figure  IV-2.  The  bandpass -f i ltered  traces  at  the 
top  of  the  figure  show  P-coda,  decaying  from  left  to  right, 
which  extend  into  and  overlap  onto  the  (presumed)  Sn  arrival 
which  arrives  approximately  15  seconds  before  the  predicted 
Sn  arrival  time.  The  problem  presented  by  these  data  is  two¬ 
fold:  (1)  without  source-to-station  calibration,  one  cannot 

be  certain  that  the  waveform  observed  at  the  predicted  arrival 
time  for  Sn  is,  indeed,  Sn;  (2)  if  the  phase  observed  is  Sn , 
there  is  a  question  as  to  where  the  phase  starts.  Applica¬ 
tion  of  the  shear  model  of  the  polarization  filter  resolves 
these  questions  as  seen  in  the  lower  set  of  traces;  the  start 
time  for  Sn  is  easily  determined,  and  further,  the  phase  can 
be  identified  with  confidence  as  a  shear  phase  on  the  basis 
of  its  being  passed  by  the  shear  model. 
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A  third  example  of  regional  waveform  extraction  is  shown 
in  Figure  IV-3.  The  bandpass -filtered  traces  at  the  top  of 
the  figure  show  shear  (Sn)  coda  extending  into,  and  mixing 
with,  waveforms  which  presumably  are  Lg,  Sg,  and  Rg .  As  with 
the  data  shown  in  the  previous  figure,  the  problem  here  is  to 
determine  whether  Lg,  Sg ,  and  Rg  can  be  identified,  and 
whether  their  start  times  can  be  determined.  The  second  set 
of  traces  in  Figure  IV-3  shows  the  result  of  applying  the 
shear  model  of  the  polarization  filter  to  the  bandpass- 
filtered  data.  The  shear  phase  Sg  can  now  easily  be  identi¬ 
fied,  and  its  start  time  is  easily  determined.  The  third  set 
of  traces  in  Figure  IV-3  shows  the  results  of  applying  the 
retrograde  elliptical  model  of  the  polarization  filter  and 
the  Love  wave  extraction  algorithm  to  the  bandpass-filtered 
data.  The  Rg  phase  can  now  be  identified  and  its  start  time 
can  be  determined  cn  the  vertical  and  radial  traces.  On  the 
transverse  trace,  Lg  is  clearly  visible. 

Figure  IV-4  shows  a  somewhat  different  extraction  prob¬ 
lem  which  must  be  resolved.  Here,  the  Pn  arrival  is  so  large 
that  is  completely  masks  later  arrivals.  Application  of  the 
compress ional  model  of  the  polarization  filter  strips  away 
the  Pn  coda,  revealing  the  Pg  arrival.  There  is  an  apparent 
compress ional  arrival  between  Pn  and  Pg  revealed  by  the 
polarization  filter.  The  origin  of  this  arrival  is  at  present 
unknown  and  should  be  investigated.  This  arrival  typifies  one 
of  the  uses  of  the  polarization  filter  —  to  extract  waveforms 
not  previously  seen  due  to  masking  coda. 

Figure  IV-5  is  presented  to  illustrate  the  effects  of 
the  polarization  filter  on  coda.  The  P  wave  shown  in  this 
figure  was  generated  by  an  event  with  an  epicentral  distance 
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FIGURE  IV- 3 

EXTRACTION  OF  Sg,  Lg ,  AND  Rg  FROM  SHEAR 
CODA  AND  SEISMIC  NOISE 
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FIGURE  1 V-4 

EXTRACTION  OF  Pg  FROM  Pn-CODA 


of  approximately  9°  and  with  a  crustal  travel  path.  The  top 
set  of  traces  shows  the  bandpass -filtered  data  while  the  sub¬ 
sequent  sets  show  the  traces  produced  by  application  of  the 
compressional,  retrograde  elliptical,  and  shear  models,  re¬ 
spectively.  After  application  of  the  compressional  model, 
the  P  coda  is  almost  completely  suppressed  while  the  maximum 
amplitudes  of  the  first  arrivals  on  the  vertical  and  radial 
components  are  reduced  less  than  five  percent.  Application 
of  the  retrograde  elliptical  and  shear  model  show  that  the  P 
coda  suppressed  by  the  compressional  model  exhibit  retrograde 
elliptical  and  shear  mode  characteristics.  This  suggests 
that  P  conversion  to  these  modes  is  taking  place  near  the 
free  surface.  Note  that  the  waveforms  on  the  retrograde  el¬ 
liptical  and  shear  model  filtered  traces  have  amplitudes  ap¬ 
proximately  one-tenth  of  the  waveforms  passed  by  the  compres¬ 
sional  model. 

Figures  IV-6  and  IV-7  illustrate  another  important  inter¬ 
pretation  of  polarization  filter  signal  extractions.  Figure 
IV-6  shows  the  P  time  gate  of  a  presumed  Honshu,  Japan,  event 
with  a  reported  epicentral  distance  of  11°;  Figure  IV-7  shows 
the  S  time  gate  for  the  same  event.  Application  of  the  com¬ 
pressional  model  of  the  polarization  filter  to  the  P  time 
gate  shows  that  the  P  arrival  actually  consists  of  an  emer¬ 
gent  Pn  arrival  (arrival  time  uncertain)  and  an  impulsive 
Pg  arrival.  Application  of  the  compressional  model  of  the 
polarization  filter  to  the  S  time  gate  shows  that  the  purport¬ 
ed  S  arrival  for  this  event  actually  consists  of  two  compres¬ 
sional  waveforms  from  a  second  event  (note  phases  tentatively 
identified  as  Pn(2)  and  Pg(2)  in  Figure  IV-7).  The  low  am¬ 
plitude  arrival  near  the  start  of  the  S  time  gate  passes  as  a 
shear  wave  (Sn  of  Figure  IV-7).  Thus,  this  reported  Honshu 
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'event'  appears,  actually,  to  be  two  events;  the  first  had 
an  epicentral  distance  of  4°-5°  (based  on  the  Pg ,  Sn  arrival 
time  difference)  while  the  second  had  an  epicentral  distance 
of  3°-4°  (based  on  the  Pn(2),  Pg(2)  arrival  time  difference). 
These  results  demonstrate  that  the  polarization  filter  can  be 
used  to  identify  phases  with  more  confidence,  and  can  result 
in  improved  event  location  and  measurement  capabilities. 

The  remaining  four  figures,  Figures  IV-8,  IV-9,  IV-10, 
and  IV-11,  illustrate  the  results  of  applying  the  polariza¬ 
tion  filter  to  mixed  signals.  Figure  IV-8  presents  the 
bandpass-filtered  vertical,  transverse,  and  radial  traces  for 
an  event  in  northeastern  China  with  an  epicentral  distance  of 
8°  (for  clarity,  the  second  half  of  the  data  window  is  pre¬ 
sented  below  the  first  half) .  To  the  left  of  the  P  tick  mark 
in  this  figure  is  seen  the  coda  of  an  earlier  event  (presumably 
consisting  of  the  Lg  and  Rg  phases).  The  first  arrival  of  the 
selected  event  can  be  seen  just  to  the  right  of  the  P  tick 
mark.  A  small  increase  in  amplitude  just  to  the  left  of  the 
Sn  tick  mark  may  indeed  be  this  phase.  Another  arrival  can 
be  seen  between  the  Sn  and  Lg  tick  marks.  Following  this  is 
the  Lg ,  Sg ,  and  Rg  train  of  waveforms.  The  amplitude  increase 
near  the  end  of  the  record  may  mark  the  arrival  of  another 
such  wavetrain.  To  resolve  and  identify  these  various  arrivals, 
the  polarization  filter  was  applied  to  the  three  time  gates 
marked  in  Figure  IV-8. 

Figure  IV-9  presents  the  results  of  applying  the  polariza¬ 
tion  filter  to  the  data  in  the  P  time  gate  of  Figure  IV-8.  The 
coda  to  the  left  of  the  signal  is  stripped  away  by  the  compres- 
sional  model,  leaving  the  clearly  defined  arrivals  Pn(l)  and 
Pg(l).  The  coda  removed  by  the  compr es sional  model  appears 
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on  the  output  from  the  retrograde  elliptical  and  shear  models, 
thereby  indicating  that  this  coda  is  the  Lg ,  Sg ,  and  Rg  wave- 
train  of  an  earlier  event. 

Figure  IV-10  presents  the  results  of  applying  the  polar¬ 
ization  filter  to  the  data  in  the  Sn-Rg  time  gate  of  Figure 
IV-8.  The  output  from  the  compress ional  model  clearly  iden¬ 
tifies  the  waveform  located  between  the  Sn  and  Lg  tick  marks 
as  a  P  wave  (identified  as  P(2)  in  Figure  IV-10).  Also  in¬ 
dicated  on  the  appropriate  traces  in  this  figure  are  the  Sn, 
Sg,  Lg,  and  Rg  arrivals  associated  with  Pn(l)  and  P  g  ( 1 )  of 
Figure  IV-9. 

Finally,  Figure  IV-11  presents  the  results  of  applying 
the  polarization  filter  to  the  data  in  the  'Post-Rg'  time 
gate  of  Figure  IV-8.  These  data  show  that  the  second  wave- 
train  of  Figure  IV-8  consists  of  the  overlapping  Lg,  Sg,  and 
Rg  waveforms  which  can  be  associated  with  arrival  P(2)  of 
Figure  IV-10.  These  arrivals,  indicated  by  Lg(2),  Sg ( 2 ) , 
and  Rg(2),  are  marked  on  the  appropriate  traces. 

In  sum,  this  section  has  illustrated  the  value  of  the 
polarization  filter  in  extracting  short-period  regional 
phases  from  seismic  noise  and  signal  coda,  in  identifying 
regional  phases,  and  in  separating  the  various  phases  of 
mixed  events. 
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SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  following  conclusions  are  drawn  from  the  work  per¬ 
formed  on  the  phase-difference  polarization  filter: 

•  Using  batch  processing,  a  broadband  bandpass  filter 
(in  this  study,  0.5-S.0  Hz)  is  more  effective  as  a 
preprocessor  to  the  polarization  filter  than  is  an 
arbitrary  narrowband  filter; 

•  Choice  of  the  processing  segment  length  used  in  the 
polarization  filter  algorithm  appears  to  have  a 
greater  effect  on  the  signal  -  to -noise  ratio  improve 
ment  of  short-period  regional  phases  than  does  the 
choice  of  the  trapezoidal  pass  window  half -width 
parameters  for  the  range  of  values  tested  (half 
widths  of  10°  to  20°)  ; 

•  The  'best'  results,  in  terms  of  signal  -  to -noise 
ratio  gain,  are  obtained  with  a  polarization  filter 
processing  segment  length  of  16  points; 

•  The  polarization  filter  provides  s ignal -to-noise 
ratio  gains  of  approximately  10  dB  for  regional 
compress ional  and  shear  waves,  and  of  approximately 
8  dB  for  Lg  and  Rg ; 

•  Direct  estimates  of  detection  capability  improve¬ 
ment  which  result  from  application  of  the  phase- 
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difference  polarization  filter  for  regional  P  and 
S  waves  were  found  to  be  approximately  0.6  m^  units 

•  In  addition  to  its  use  as  a  signal  extractor,  the 
phase-difference  polarization  filter  can  be  used 
to  separate  mixed  signals  and  to  identify  phases. 


B.  RECOMMENDATIONS  FOR  FUTURE  USE 

The  following  recommendations  should  be  considered  in 
any  future  work  involving  the  phase -dif f erence  polarization 
filter : 

•  Implement  the  phase-difference  polarization  filter 
in  an  interactive  processing  environment; 

•  Cascade  the  phase-difference  polarization  ^ilter 
with  the  previously-developed  particle-motion 
polarization  filter  (Strauss,  1978)  which  tests  on 
the  angle  of  incidence  to  create  a  total  model  of 
bodywave  polarization; 

•  Test  the  phase-difference  polarization  filter  on 
data  recorded  at  a  different  station  in  order  to 
check  for  regional  dependence  of  program  parameters 
and  in  order  to  estimate  better  the  Lg  and  Rg  ex¬ 
traction  results. 

NOTE: 

Advanced  applications  of  polarization  filtering  should 
address  a  more  complete  interpretation  of  phase  arrival  times 
and  of  the  relation  of  these  arrivals  to  local  crustal  struc¬ 
ture.  Also,  the  relative  excitation  and  frequency  of  phases 
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can.  be  used  to  determine  information  on  source  depth  and 
source  mechanism.  This  important  application  can  be  facil¬ 
itated  by  further  development  of  interactive  waveform  extrac¬ 
tion  combined  with  synthetic  studies  which  would  be  used  to 
predict  the  observed  wavetrain  associated  with  specific  source 
region-to-station  paths.  Results  shown  here,  with  a  few  of 
the  better  known  phases,  indicate  interesting  possibilities 
for  applying  this  approach  to  precise  seismic  interpretation 
of  observed  phases. 

Criteria  for,  and  the  precision  of,  source  direction 
and  distance  determinations  should  be  thoroughly  evaluated. 

The  combined  use  of  the  phase  associations  with  source  dis¬ 
tance,  direction,  and  depth  estimates  may  provide  the  only 
feasible  basis  for  focal  determinations  and  discrimination 
of  small  regional  events. 
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APPENDIX  A 

THE  POLARIZATION  FILTER  PROGRAM 


This  appendix  describes  the  operation  of  the  polariza¬ 
tion  filter  program  as  it  is  presently  constituted. 

Figure  A-l,  a  flow  diagram  of  the  polarization  filter 
program,  shows  the  steps  taken  in  processing  data.  Data  to 
be  processed  consist  of  a  maximum  of  409.6  seconds  of  short- 
period  data  (4096  data  points  sampled  at  10  points  per  sec¬ 
ond),  the  horizontal  components  of  which  have  been  rotated 
from  their  recorded  north-south,  east-west  orientation  to 
a  transverse,  radial  orientation. 

Once  entered  into  the  program,  the  data  are  first  band¬ 
pass  filtered  to  provide  a  standard  against  which  to  compare 
the  polarization  filter  output.  Next,  the  RMS  noise,  peak 
signal  amplitudes,  and  signal- to-noise  ratios  for  the  com- 
pressional,  shear,  and  surface  waves  are  computed  from  the 
bandpass  filtered  data.  The  bandpass  filtered  data  next 
serve  as  input  to  the  polarization  filter  algorithm  in  order 
to  achieve  maximum  SNR  gain  (this  filter  cascading  approach 
was  developed  by  Lane  (1977)  and  Strauss  (1978)). 

The  next  step  in  the  processing  sequence  is  to  check  a 
set  of  flags  input  to  the  program  to  determine  the  type  of 
waveform  (compressional ,  retrograde  elliptical,  shear,  or 
prograde  elliptical)  the  polarization  filter  algorithm  is 
to  pass.  Each  of  the  four  flags  either  has  a  value  of  zero 
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FIGURE  A- 1 

GENERALIZED  FLOW  OF  PHASE- DIFFERENCE  POLARIZATION  FILTER 


denoting  no  polarization  filtering  is  to  be  performed  for 
this  model,  or  one  denoting  polarization  filtering  is  to  be 
performed.  Thus,  if  the  first  flag  is  set  to  one,  the  com- 
pressional  waveform  model  of  the  polarization  filter  is  ap¬ 
plied.  Upon  completion  of  waveform  processing,  the  second 
flag  is  tested  to  determine  whether  the  retrograde  elliptical 
waveform  model  is  to  be  applied.  Following  this,  the  third 
flag  is  tested  to  determine  whether  the  shear  waveform  model 
is  to  be  applied,  and,  finally,  the  fourth  flag  is  tested  to 
determine  whether  the  prograde  elliptical  waveform  is  to  be 
applied  . 

Once  the  program  has  determined  that  a  particular 
model  is  to  be  applied,  the  controlling  parameters  of  the 
polarization  filter  algorithm  (i.e.,  the  number  of  points 
per  processing  segment,  and  the  width  of  the  pass  window) 
are  set  to  their  pre -determined  values.  The  bandpass  filter¬ 
ed  data  are  then  polarization  filtered  according  to  the 
selected  model.  From  the  filtered  data,  the  RMS  noise,  peak 
signal  amplitude,  and  SNR  are  computed,  and  the  filtered  data 
are  written  on  a  plot  tape.  The  program  then  restores  the 
bandpass  filtered  data  and  loops  back  for  the  next  waveform 
model.  The  RMS  noise,  peak  signal  amplitude,  and  signal-to- 
noise  ratio  computed  before  and  after  polarization  filtering 
are  printed  out  to  allow  computation  of  noise  suppression, 
signal  degradation,  and  signal -to-noise  ratio  improvement 
due  to  the  polarization  filter. 

Figure  A- 2  is  a  detailed  flow  diagram  of  the  polariza¬ 
tion  filter  algorithm.  To  begin,  the  first  data  segments 
(say,  32  points  of  the  vertical,  radial,  and  transverse 
traces)  are  moved  into  a  processing  array.  Next,  the  mean 


ENSCO,  INC. 


A  -  3 


FIGURE  A- 2 

DETAILED  FLOW  OF  POLARIZATION  FILTER  ALGORITHM 


ENSCO,  INC. 


A  -  4 


t 


1 


of  each  data  segment  is  removed,  the  processing  array  is  pad¬ 
ded  with  zeros  to  eight  times  the  segment  length,  and  the  re¬ 
sulting  zero-mean  data  are  Fourier  transformed.  The  algorithm 
then  enters  a  loop  on  the  number  of  discrete  frequencies.  For 
each  frequency,  the  difference  in  phase  between  the  radial  and 
vertical  components  is  computed.  This  phase  difference  is 
then  compared  to  the  trapezoidal  pass  window  (shown  in  Figure 
A- 3) ,  which  is  centered  on  the  expected  phase  difference  for 
the  waveform  model  under  consideration.  The  filter  weight  for 
this  frequency,  as  determined  by  this  comparison,  is  then  ap¬ 
plied  to  the  data  at  that  frequency.  Specifically,  it  is  ap¬ 
plied  to  three  traces  for  the  compress ional  and  shear  models, 
and  to  the  vertical  and  radial  components  for  the  retrograde 
and  prograde  elliptical  models. 

The  purpose  in  applying  the  filter  weights  to  the  trans¬ 
verse  component  in  the  compress ional  case  is  to  pass  any 
compress ional  energy  present.  Significant  amounts  of  cora- 
pressional  energy  on  this  component  indicate  an  off-azimuth 
arrival;  i.e.,  the  waveform  being  extracted  is  not  from  the 
expected  source.  Applying  the  filter  weights  to  the  trans¬ 
verse  component  in  the  shear  case  permits  extraction  of  S^, 
which  arrives  at  the  same  time  as  Sy.  For  the  retrograde 
and  prograde  elliptical  models,  the  transverse  data  segment 
is  processed  by  a  separate  algorithm  which  is  designed  to 
extract  Love  waves.  This  algorithm  tracks  the  incoming  Love 
wave  in  azimuth  and  rejects  all  energy  with  approach  azimuths 
outside  a  preset  range  centered  on  the  great  circle  azimuth 
(see  Strauss ,  19  76) . 

After  completion  of  the  loop  on  frequencies ,  the 
polarizat ion- filtered  segments  are  inverse  fast  Fourier 


ENSCO,  INC. 


A  -  5 


transformed  and  saved.  Ensuing  segments  each  overlap  the 
preceding  segment  by  50  percent.  After  processing,  the  ove 
lapping  data  points  are  scaled  and  summed  to  construct  the 
polarization- filtered  time  traces. 
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